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Though the notion of phase synchronization has been well studied in chaotic dynamical systems 
without delay, it has not been realized yet in chaotic time-delay systems exhibiting non-phase coher- 
ent hyperchaotic attractors. In this article we report the first identification of phase synchronization 
in coupled time-delay systems exhibiting hyperchaotic attractor. We show that there is a transition 
from non-synchronized behavior to phase and then to generalized synchronization as a function 
of coupling strength. These transitions are characterized by recurrence quantification analysis, by 
phase differences based on a new transformation of the attractors and also by the changes in the Lya- 
punov exponents. We have found these transitions in coupled piece-wise linear and in Mackey-Glass 
time-delay systems. 
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Synchronization is a natural phenomenon that one en- 
counters in daily life. Since the identification of chaotic 
synchronization Q, |M |1{ , several papers have appeared 
in identifying and demonstrating basic kinds of synchro- 
nization both theoretically and experimentally (cf . 0, ) . 
Among them, chaotic phase synchronization (CPS) refers 
to the coincidence of characteristic time scales of the 
coupled systems, while their amplitudes of oscillations 
remain chaotic and often uncorrelated. Phase synchro- 
nization (PS) plays a crucial role in understanding the 
behavior of a large class of weakly interacting dynam- 
ical systems in diverse natural systems. Examples in- 
clude circadian rhythm, cardio-respiratory systems, neu- 
ral oscillators, population dynamics, electrical circuits, 
etc 



The notion of CPS has been investigated so far in os- 
cillators driven by external periodic force 0, II] , chaotic 
oscillators with different natural frequencies and/or with 
parameter mismatches 0, [HI |3 ) arrays of coupled 
chaotic oscillators llJ a nd also in essentially differ- 
ent chaotic systems |l5l Il6l|. On the other hand PS in 
nonlinear time-delay systems, which form an important 
class of dynamical systems, have not yet been identified 
and addressed. A main problem here is to define even 
the notion of phase in time-delay systems due to the in- 
trinsic multiple characteristic time scales in these sys- 
tems. Studying PS in such chaotic time-delay systems 
is of considerable importance in many fields, as in un- 
derstanding the behavior of nerve cells (neuroscience), 
where memory effects play a prominent role, in patho- 
logical and physiological studies, in ecology, in lasers etc 



iiBEIIllTiailil. 

In this paper, we report the first identification of phase 
synchronization (PS) in nonidentical time-delay systems 
in the hyperchaotic regime with non-phase coherent at- 
tractors with unidirectional nonlinear coupling. We will 
show the entrainment of phases of a coupled piecewise 
linear time-delay system for weak coupling from the 
non-synchronized state. Phase is calculated using the 
Poincare method 0, 0| after a new transformation of at- 
tractors of the time-delay systems, which looks then like 
a smeared limit cycle. The existence of PS and general- 
ized synchronization (GS) in coupled time-delay systems 
is characterized by recently proposed methods based on 
recurrence quantification analysis and also in terms of 
Lyapunov exponents of the coupled time-delay systems. 

We first consider the following unidirectionally coupled 
drive xi{t) and response ^2 (^) systems, which we have 
recently studied in detail in [23, llj, Hj] , 

xi{t) = -axi{t) + bif{xi{t-T)), (la) 

X2{t) = -aX2it)+b2fix2{t-T))+b3f{xi{t-TX)\h) 

where 61,62 and 63 are constants, a > 0, t is the delay 
time and f{x) is the piece- wise linear equation of the 
form 



/(^) = 




X < -4/3 
-4/3 <x< -0.: 
-0.8 < a; < 0.8 
0.8<x < 4/3 

X > 4/3. 



(2) 
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We have chosen the values of parameters as a = 1.0, 61 — 
1.2,62 = 1-1 and T — 15, which are outside the region 
of complete, lag and anticipatory synchronizations dis- 
cussed in |23,|23|- For this parametric choice, in the ab- 
sence of coupling, the drive Xi{t) and the response X2{t) 
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FIG. 1: Phase synchronization between the systems and 
(03). (a) The non- phase coherent hyperchaotic attractor of 
the uncoupled drive 1^). (b) Transformed attractor in the 
xi{t + r) and z{t + r) space along with the Poincare points 
represented as open circles, (c) The phases of the drive system 
{4>i), the response system (01) calculated from the new state 
variable z{t + r) and their difference (A^) for 63 = 1.5. 



systems evolve independently. Further in this case, the 
drive xi{t) exhibits a hyperchaotic attractor (Fig. ^) 
with five positive Lyapunov exponents (See ||2^ for the 
spectrum of Lyapunov exponents) and the response X2 (t) 
has four positive Lyapunov exponents, i.e. both subsys- 
tems are qualitatively different (due to 61 7^ 62)- The 
parameter 63 is the coupling strength of the unidirec- 
tional nonlinear coupling ^p), while the parameters 61 
and 62 play the role of parameter mismatch resulting in 
nonidentical coupled time-delay systems. 

Now the important questions we encounter are whether 
PS exists in the time-delay system ^ when the coupling 
is included (63 > 0) and, if so, how to characterize the 
possible transition to PS in such systems which possess 
in general highly non-phase coherent attractors having 
a broad-band power spectrum. In low-dimensional sys- 
tems, a few methods 0,0 have been developed to define 
phase in phase coherent chaotic attractors, which have 
a dominant peak in the power spectrum. Definition of 
phase is not so clear in non-coherent chaotic attractors, 
in particular in high-dimensional systems having broad- 
band power spectra such as time-delay systems. Meth- 
ods to calculate phase of non-coherent attractors of time- 
delay systems is not readily available. One approach to 
calculate the phase of system is based on the concept 
of curvature |25j |. which is often used in low-dimensional 
systems. However, we find that this procedure does not 
work in the case of time-delay systems in general, and 
in particular for Eqs. We present here three other 
approaches to study PS in systems as JQ): 

i) We introduce a new transformation to successfully 
capture the phase in the present problem. It transforms 
the non-phase coherent attractor (Fig.^) into a smeared 
limit cycle-like form with well defined rotations around 



one center (Fig. ^). This transformation is performed 
by introducing the new state variable 

Z{t + T) = Xl{t)xi{t + T)/Xl{t + T), (3) 

where f is the optimal value of delay time to be chosen 
(so as to rescale the original non-phase coherent attractor 
into a smeared limit cycle-like form), we plot the above 
attractor (Fig.^) in the {xi{t + T), z{t + T)) phase space. 
The functional form of this transformation has been iden- 
tified by generalizing the transformation used in the case 
of chaotic atractors in the Lorenz system We find 
the optimal value of f to be 1.65. The above transfor- 
mation is obtained through a suitable functional form 
(along with a delay time f), so as to unfold the origi- 
nal attractor (Fig. Q^) into a phase coherent attractor. 
Now the attractor (Fig. ^p) looks indeed like a smeared 
limit cycle with nearly well defined rotations around a 
fixed center. Hence, we can calculate the phase using 
the Poincare method |3, from the rescaled attractor. 
The Poincare points are shown as open circles in Fig.^. 
The phases of both the drive (t) and the response ^| (t) 
systems calculated from the new state variable z(t + t) 
are shown in Fig. ^ along with their phase difference 
A(j) = 0f (i) — 4>2{t) for the value of the coupling strength 
63 = 1.5, showing a high quality PS. This strong bound- 
edness of the phase difference obtains for 63 > 1.382. 
(Note that the transformed attractor (Fig. ^) docs not 
have any closed loops as in the case of the original at- 
tractor (Fig. ^). If it is so, such closed loops will lead 
to phase mismatch, and one cannot obtain exact match- 
ing of phases of both the drive and response systems as 
shown in Fig. ^). 

ii) Next, we analyze the complex synchronization phe- 
nomena in the coupled time-delay systems JQl by means 
of ver y re cently proposed methods based on recurrence 
plots [2g. These methods help to identify and quan- 
tify PS (particularly in non-phase coherent attractors) 
and GS. For this purp ose, the generalized autocorrela- 
tion function P{t) |2g was introduced 

N-t 

^w^]^E©(^-ii^'-^»+*ii)' (4) 

i=l 

where 6 is the Heaviside function, Xi is the ith data 
of the system X, e is a predefined threshold. ||.|| is the 
Euclidean norm and N is the number of data points. 
Looking at the coincidence of the positions of the maxima 
of P{t) for both systems, one can qualitatively identify 
PS. 

A criterion to quantify PS is the cross correlation coef- 
ficient between the drive, Pi{t), and the response, P2{t), 
which can be defined as Correlation of Probability of Re- 
currence (CPR) 

CPR^ {Pi{t)P2{t))/aia2, (5) 

where Pi_2 means that the mean value has been sub- 
tracted and ai^2 are the standard deviations of Pi{t) and 
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FIG. 2: Generalized autocorrelation functions of both the 
drive P\{t) and the response P2{t) systems, (a) Non-phase 
synchronization for 63 = 0.6, (b) Phase synchronization for 
63 = 1.5 and (c) Generalized synchronization for 63 = 2.3. 



P2{t) respectively. CPR « 1 indicates that the systems 
are in complete PS, whereas for non-PS one obtains low 
values of CPR. 

To characterize GS, the authors of '26l| proposed the 
first index as the Joint Probability of Recurrences (JPR), 



JPR = 



1 

JV2 



\X,-X,\\)Q{ey 



\Y,~YA\)- 



1 - RR 



(6) 



where RR is rate of recurrence, Cx and Cy are thresholds 
corresponding to the drive and response systems respec- 
tively. RR measures the density of recurrence points and 
it is fixed as 0.02 [26|. JPR is close to 1 for systems in 
GS and is small when they are not in GS. The second 
index depends on the coincidence of probability of re- 
currence, which is defined as Similarity of Probability of 
Recurrence (SPR), 



SPR - 1 - {{Pi{t) - P2{t)f)/aia2. 



(7) 



SPR is of order 1 if both the systems are in GS and 
approximately zero if they evolve independently. 

Now, we will apply this concept to the original (non- 
transformed) attractor (Fig. ^). We estimate these re- 



currence based measures from 5000 data points after suf- 
ficient transients with the integration step h = 0.01 and 
sampling rate At = 100. The generalized autocorre- 
lation functions Pi{t) and P2{t) (Fig. for the cou- 
pling 63 = 0.6 show that the maxima of both systems 
do not occur simultaneously and there exists a drift be- 
tween them, so there is no synchronization at all. This 
is also reflected in the rather low value of CPR = 0.381. 
For 63 G (0.91,1.381), we observe the first substantial 
increase of recurrence reaching CPR w 0.5 — 0.6. Look- 
ing to the details of the generalized correlation functions 
P(t), we find that now the main oscillatory dynamics 
becomes locked, i.e. the main maxima of Pi and P2 co- 
incide. For 63 e (1.382,2.2) CPR reaches almost 1, i.e. 
now all maxima of Pi and P2 are in agreement and this is 
in accordance with strongly bounded phase differences. 
This is a strong indication for PS. Note, however that the 
heights of the peaks are clearly different (Fig. I^b)- The 
differences in the peak heights exhibit that there is no 
strong interrelation in the amplitudes. Further increase 
of the coupling (here 63 = 2.3) leads to the coincidence of 
both the positions and the heights of the peaks (Fig. [3:) 
referring to GS in systems This is also confirmed 
from the maximal values of the indices JPR — 1 and 
SPR = 1, which is due to the strong correlation in the 
amplitudes of both systems. The transition from non- 
synchronized to PS and then to GS is characterized by 
the maximal values of CPR, SPR and JPR (Fig.OjD). As 
expected from the construction of these functions, CPR 
refers mainly to the onset of PS, whereas JPR quanti- 
fies clearly the onset of GS. The existence of GS is also 
confirmed using the auxiliary system approach ^2^. 

iii) The transition from non-synchronization to PS is 
also characterized by changes in the Lyapunov exponents 
of the coupled time-delay systems The spectrum of 
the nine largest Lyapunov exponents of the coupled sys- 
tems is shown in Fig. from which one can find that 
the Lyapunov exponents corresponding to the response 
system become negative from the value of the coupling 

strength 63 > 0.9, where the transition to PS occurs, 

(2) 

RRexcept for the largest Lyapunov exponent Xmax which 
continues to remain positive. This is a strong indication 
that in this rather complex attractor the amplitudes be- 
come somewhat interrelated already at the transition to 
PS (as in the funnel attractor [2^). It is interesting to 

note that the Lyapunov exponents of the response system 
(2) 

Xi (other than Xmax) are changing already at the early 
stage of PS (63 « 0.91), where the complete PS is not yet 
reached. This has been also observed for the onset of PS 
in phase-coherent oscillators IJj. 

We have obtained the same results for different sam- 
pling intervals At and for various values of delay time 
T. We have also identified this transition to PS and to 
GS in the coupled Mackey-Glass systems 0, [23, and 
also in for linear coupling (unidirectional) using these 
above three approaches (these results will be published 
in a more comprehensive paper). 

In conclusion, we have identified for the first time the 
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FIG. 3: (a) Spectrum of first nine largest Lyapunov expo- 
nents of the coupled systems (0 as a function of coupling 
strengtii 63 € (0, 3) and (b) Spectrum of CPR and JPR as a 
function of coupling strength bg G (0, 3). 



existence of PS in coupled time-delay systems in the hy- 
perchaotic regime with highly non-phase coherent attrac- 
tors. We have shown that there is a typical transition 
from non-synchronized state to PS for weak coupling and 
in the range of strong coupling there is a transition to GS 
from PS. We have also identified a suitable transforma- 
tion, which works equally well for Mackey-Glass system 
(having a more complex hyperchaotic attractor), to cap- 
ture the phase of the underling attractor. We have also 
characterized the existence of PS and GS in terms of 
recurrence based indices like generalized autocorrelation 
function P{t), CPR, JPR and SPR and quantified the 
different synchronization regimes in terms of them. The 
above transition is also confirmed by the changes in the 
Lyapunov exponents. We have pointed out the existence 
of PS in coupled Mackey-Glass systems as well. The re- 
currence based technique as well as the new transforma- 
tion are also appropriate for the analysis of experimental 
data, i.e. we expect experimental verification of these 
findings. 
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